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Abstract. It is well-known that the functions Vm (x) = r , /n°° " du 

arise naturally in the study of atoms in strong magnetic fields where they can 
be regarded as one-dimensional regularizations of the Coulomb potential. For 
many-electron atoms consideration of the Pauli principle requires convex com- 
binations of such potentials and interactions of the form V m ( ^^g^ ) • We 
summarize the results of a comprehensive study of these functions including 
recursion relations, tight bounds, convexity properties, and connections with 
confluent hypergeometric functions. We also report briefly on their application 
in one-dimensional models of many-electrons atoms in strong magnetic fields. 



1. Introduction 

It is well-known that systems in strong magnetic fields behave like systems in 
one-dimension, i.e., a strong magnetic field confines the particles to Landau orbits 
orthogonal to the field, leaving only their behavior in the direction of the field 
subject to significant influence by a static potential. Motivated by this general 



principle and the work of Lieb, Solovej and Yngvason [LSY| on atoms in extremely 



strong magnetic fields, Brummelhuis and Ruskai [BR1] initiated a study of models 
of atoms in homogeneous strong magnetic fields in which the 3-dimensional wave- 
function has the form 

(1-1) *(ri,r2 . . .r„) = ip(xi . . .x n ) T(j/i, z x ,y 2 , z 2 , ■ ..y„, z n ) 

where T lies in the projection onto the lowest Landau band for an N-electron system. 
We follow the somewhat non-standard convention of choosing the magnetic field in 
the x-direction, i.e, B = (B. 0, 0) where B is a constant denoting the fields strength, 
in order to avoid notational confusion with the nuclear charge Z . 
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Such models lead naturally to one dimensional regularizations of the Coulomb 
potentials of the form 



V m (x) 
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Recognition that such potentials arc impor- 



(1.2) 
(1.3) 
(1.4) 

where 7m(r, #) 

tant goes back at least to Schiff and Snyder SSI in 1939 and played an important 
role in the Avron, Herbst and Simon study [AHS| of hydrogen. Recently, Ruskai 
and Werner [RW] undertook a detailed study of these potentials, proving the im- 
portant property of convexity of 1/V m as well as a number of other useful properties. 
The primary purpose of this note is to give a summary of these results in Sections 
[| and [| Before doing that, we briefly discuss one-dimensional models of atoms in 
strong magnetic fields in Section |^ and their implications for the maximum negative 
ionization problem in Section |3[ 



(2.1) 



2. Atoms in Strong Magnetic Fields 

The Hamiltonian for an N electron atom in a magnetic field B is 

1 



H(N, Z,B) 



N 

E 



Z 



E 



where A is a vector potential such that V x A = B. The ground-state energy of 
H(N, Z, B) is given by 



(2.2) 



Eo(N, Z,B) 



inf (H(N,Z,BW,y) 

ll*ll=i 



Let En° ni (N, Z, B) denote the corresponding minimum restricted to functions of the 



form (1.1). For extremely strong fields, it was shown in [LSY| that Eq/E\ 



473 



T^conf 
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with N/Z fixed. 



as B/Z 

In | BRl |and BR2 | we consider two special cases of (LI). We write the Landau 
state with angular momentum m in the x-direction in the form 



(2.3) 



B (m+l)/2 
V7T 



■c 
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where C, = y + iz. Then our two special cases can be described as follows: 

Zero model: In this case, we make the extremely simple assumption that T 

is simply a product of Landau states with m = 0, i.e., T = Jlfcli lo{Vk, 
Slater model: In this case we assume that T is an antisymmetrized product 

of Landau states with m = 0, 1, . . . N— 1, i.e., T = -^=[y A 71 A . . . A^n-i]- 

Although the first model is somewhat unrealistic, its simplicity makes it amenable 
to detailed analysis which yields insight into the general situation. The second 
model corresponds to the physically reasonable assumption that T is a Slater de- 
terminant. In this case, the required antisymmetry of the wave function is inherent 
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in our assumptions on T and the one dimensional function ip is symmetric, i.e., the 
electrons behave like bosons in one dimension. 
It is straightforward to show that 



(2.4) 



E£ onl (N,Z,B) = VB 



inf {h(N,Z,M)i/),ip) 

•0||=i 



NB 



where have scaled out the field strength B so that 

1 d 2 



(2.5) 



h(N, Z, M) = 



N 

E 



M dx 2 v 3 ' 



j<k 

and the only remnant of the magnetic field is in the "mass" M = B~ x l 2 and the 
effective one-dimensional potentials V and W will be defined below for each model. 



" " T2 V ° 



For the zero model, one easily finds BRl| that 

(2.6) V(x) = V (x) and W(u 
Note Vo(x) « t~t for large |jc|. Thus, for large separations 

(2.7) W(u-v) 



\u — v\ 



1 \/2 



f 



V2\ 



However, not only is the singularity removed at u = v, W(0) is smaller by a factor 
of 4g than V(0). This means that if two electrons are simultaneously near the 
nucleus, the price paid from the electron-electron repulsion is smaller than that 
gained from the electron-nuclear attraction. Although this effect seems to play an 
important role in binding additional electrons, it may be partially offset by the 
price paid in kineti c energy if one attempts to constrain both electrons near the 
nucleus. Sec [BRll for further discussion. 



For the Slater model it can be shown [BR2] that 



N-l 



(2.8) 
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V(x) = V^(x) = ^J2 V ^ and 
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W(u-v) = —= 
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where Cj > V j and J2j c j 
bination of V m with odd m 



1 so that the effective interaction is a convex com- 
1,3,... 2N — 1, albeit with the same scaling as 



in (2.7). Note that the convex sum in (|2.9| ) above includes contributions from V m 
with m > N. Properties (b) and (d) of Section [| imply that V m (Q) is decreasing 
in rn. Therefore, one expects a decrease in the electron-electron repulsion W in 
addition to that from the factor of However, delicate combinatorics would be 
needed to verify this exactly. 

Some obvious variations on these models are possible and discussed briefly in 

|bri| , 



BR2]. It is interesting to note that if T 



T-fN 



{Vk, Zk) with to odd, 

then the convex sum analogous to ( |2.9| ) contains only terms V^j with even subscript. 
It is also worth noting that 



(2.10) 



lim 

/3^oo log p 



V m {fix) = 5{x) 
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in the sense of tempered distributions. This implies that the potentials V and W 
which occur in our models have an analogous delta potential behavior as f3 — > oo. 
The proof ||d| of ( pTTo| ) uses the Fourier transform (property (k) in Section |]) 



of the potentials V m (x), particularly the observation that V m ((,) has a logarithmic 
singularity at £ = 0. Similar limiting behavior (with /3 = y/B) was observed by 
LSY] for potentials in the three-dimensional Hamiltonian (2.1). It can also be 
shown [ BD[ that if ( |2.5| ) is appropriately rescaled and the potentials replaced by 
the corresponding delta potentials, the result is a one-dimensional Hamiltonian 
whose semi-classical limit (on bosonic wave functions) as Z — > oo is precisely that 



given by the (fermionic) hyperstrong functional in LSY]. Since their functional was 
shown to describe the Z, B/Z 3 — ► oo limit of the three-dimensional Hamiltonian 



(2.1), this provides additional justification for our models. Even the simple-minded 
zero model has the correct asymptotic behavior. 

3. Maximum Negative Ionization 

In the absence of a magnetic field, one expects that the maximum number 
of electrons a nucleus with charge Z can bind is N max (Z) = Z + 1 or Z + 2. 
However, only the somewhat weaker result of asymptotic neutrality has been proved 
rigorously [LSST]. If electrons behave like bosons, asymptotic neutrality does not 
hold and iV max behaves asymptotically roughly like \.2\Z. (See Sol | for details and 
references to earlier work on bosonic atoms.) In [Lb] Lieb gave a simple argument 
which showed that A max < 2Z + 1, independent of particle statistics. Thus, it may 



seem somewhat surprising that [LSY] showed that for atoms in extremely strong 
magnetic fields 

A max (Z) 



(3.1) 



lim inf 

Z, B/Z 3 ^c 



> 2. 



The study of one-dimensional models in [BR1 was initiated, in part, by the hope 
of proving an asymptotic upper bound of the form iV max < 2Z as B,Z — > oo. 
Although we did not succeed in proving such a bound, even for our simplified one- 
dimensional models, we believe that they offer considerable insight into both the 
reasons for binding an "extra" Z electrons and the reasons why the localization 
techniques developed to bound iV max fail in the strong field case. 

It is generally believed that enhanced binding occurs in strong magnetic fields 
because the field confines the electrons in two dimensions and effectively reduces 
the atom to a one-dimensional system. Although there is some truth to this, it was 
shown in [LSY] that atoms do not become truly one-dimensional unless B > Z 3 
and the field strength is greater than anything seen on earth. (Sufficiently strong 
magnetic fields do exist on the surface of neutron stars, making this analysis of some 
interest in astrophysics.) Moreover, the binding enhancement achieved by making 
the system effectively one-dimensional can only account for small effects, such as 



the fact [AHS] that singly negative ions always have infinitely many bound states 



in a magnetic field. It cannot account for the binding of an additional Z electrons. 



The results in [BR1 suggest that the primary mechanism for binding ad- 
ditional electrons in strong fields is the fact that the effective reduction in the 
strength of the electron-electron repulsion permits two electrons to be near the 
nucleus simultaneously. However, the one-dimensional confinement also delocalizes 
the electron. This effect is seen in the Hamiltonian h{N,Z,B^ x / 2 ) given by ( |2.5P 
where the effective mass is M = B^ 1 ^ 2 so that in strong fields the electrons behave 
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like extremely light particles. The uncertainty principle then implies that trial wave 
functions which localize the electrons cannot yield bound states. 



Since Lieb's strategy [Lb| for finding an upper bound on -/V max (Z) does not 
require an explicit localization, it might seem well-suited to atoms in strong mag- 
netic fields. However, Lieb's method actually has an implicit localization (which is 



based on an idea of Benguria [Ben] for spherically symmetric atoms) for which the 



localization error is zero in three dimensions. However, as explained in |BR1|, the 
localization error is necessarily non-zero in one-dimension. (This is a consequence 
of the fact that non-positive potentials always have at least one bound state in one 
dimension. Thus, the phenomenon of enhanced binding in one dimension actually 
contributes to the derealization of the electrons.) Using Lieb's method for the zero 



model, we were only able to show in | BR1 that iV max (Z, B) < 2Z + 1 + c\/~B for 
an explicit constant c. In the interesting case B = 0(Z 3 ), this yields a bound of 
the form JV max < 2Z + cZ 3 / 2 , rather than a linear one. 

Surprisingl y, one can get a better bound using the Ruskai-Sigal localization 



method. (See Rusk | for a summary.) For both the zero model and the Slater 



model, we can prove the following result. 

Theorem 3.1. Let N max (Z, B) be the maximum number of electrons for which 



the Hamiltonian (2.1) has a bound state, and assume that the potentials V and W 



have either the form (2.6) corresponding to the zero model or the form (2.t) and 
( 2. i ) corresponding to the Slater model. Then for every a > and (3 > there is 
a constant C a p such that 

(3.2) N max {Z,B) < C aP Z 1+a BP 

where a, /3 can be arbitrarily small and, in the case of the Slater model B > Z 3+1 
for some 7 > 0. 

This result can be improved slightly to 

(3.3) N max (Z,B) < [Z(\ogZf + Z\ogZ(\ogB) 1+ "] 

where, as above, to > can be arbitrarily small and in the case of the Slater model 
B > Z 3+1 . Because the electrons in the one-dimensional model are essentially 
bosonic, this is the best that one can hope to achieve with the Ruskai-Sigal method. 

In the case of the Slater model, the Landau level portion of the wave function 
T is antisymmetric. Hence, the one-dimensional part of the wave function ip must 
be symmetric. Even for the zero model, it is physically reasonable to treat the 
electrons in the one-dimensional model as essentially bosonic. In the Ruskai-Sigal 
method, the system is divided into a small "inner" ball in which binding is precluded 
because the electrons are confined to a small region, and an "outer" ball in which 
the localization error becomes negligible as Z — > 00. For bosonic systems, one 
can always squeeze the electrons closer together, yielding a smaller cut-off p than 
for fermions. This feature is the only factor which precludes extending the proof 



of asymptotic neutrality in [LSST| to bosonic atoms. This demonstrates that the 
localization error is not simply a technical artifact, but a reflection of a real physical 
effect. 

For atoms in strong magnetic fields, the cut-off radius p is not small. Instead 
p ~ N VBZ- 2 (log ^)~ 2 which grows with B. For B = 0(Z 3 ) and N = O(Z), 
roughly speaking (i.e., ignoring the log term) p ~ Z 1 / 2 ~ B 1 ^ 6 . Thus, the lo- 
calization method can be used to obtain a (non-optimal) upper bound on N max 
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despite the fact that the electrons are highly delocalized and the size of the "inner" 
region becomes arbitrarily large as B — > oo. However, the non-optimal bounds 
above are probably the best one can expect from configuration space localization. 
It seems likely that a proof of better upper bounds will require the use of phase 
space localization techniques. 

4. Properties of V m (x) 

The functions V m (x) are even and well-defined for x > 0. Although the primary 
interest in physical applications is for integer m > 0, it is easy to see from the form 



(1.2) that they can be extended to complex m with 5R(m) > — 1. For 5ft(m) > — \ 
they are also well-defined for x = 0. In this note we restrict ourselves to non- 
negative x and real m > — 1. 

It is convenient to define V-i(x) — j^r and note that this is justified in the 

sense that lim x~ 1 V m (x) — 1 for all x > 0. 
m—-i+ 

We now summarize the properties of V m (x) for x £ (0, oo). Unless otherwise 



stated to is real and to > — 1. For proofs and further discussion, sec [KW|. 
Summary of Properties of V m (x): 

a) V m (x) satisfies the inequality = > V m (x) > —^^=^^= where the 

Vi 2 + m y x 2 + m + 1 

upper bound holds for m > and the lower for to > — 1. 

b) V m (x) is decreasing in to. In particular, V m +i(x) < V m (x) < —. 

x 

c) The expression mV m (x) is increasing in to > — 1. 

d) For to > —1/2, the definition of V m (x) can be extended to x — and 

, r(m+|) 
V^(0) = 



T(to+ 1) ' 

For integer m, this becomes 

v - ( 2w ) ! ^_ 1-3-5...(2to-1) 
22m(m!)2 VT- 2 -4- 6... (2m) 

while for large to Stirling's formula implies 

'm-i\ m /e\i/2 1 



\ m J \m/ 



which is consistent with property (a). 

e) For all to > 0, V rn satisfies the differential equation 

f) For each fixed to > 0, V m (x) is decreasing in x. 

g) For a > 0, the expression aV m (ax) increases with a. Hence aV m (ax) > V(x) 
when a > 1 and aV m (ax) < V{x) when a < 1. 

h) Vq(x) is convex in x > 0; however, V^(a;) is not convex when to > i. 

i) For integer to, 1/V^n(a;) is convex in x > 0. 

j) For integer to, the ratio V m +\(x) /V m {x) is increasing in x > 0. 

k) The Fourier transform is given by 

^ i poo Am+l poo rn — s 

s ^ L w e "" S ^ - VST L WT^ ds 
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1) For large x, it follows from property (a) that 

to 1 . . m + 1 

< - - V m (x) < 



2(x 2 + to) 3 / 2 a; 



2a; 3 



while (1.3) yields the asymptotic expansion 

1 to + 1 

a; 2a; 3 



V m (x) 



3(to + 2)(to+ 1) ^ / 1 



The lower bound in (a) was proved earlier (at least for integer to) by Avron, 
Herbst and Simon [ AHS[ . Properties (b) and (c) imply that V m (x) is decreasing 
in m, while mV m (x) is increasing; this gives an indication of the delicate behavior 
of V m . The differential equation (e) can be verified using integration by parts in 
(|Q|). Property (f) follows directly from (b) and (e). Property (g) follows from 

( |l.3| ) and the observation that — is increasing in a. ft is useful in analyz- 

V a 2 x 2 + u 

ing W(u — v) since it implies (^j^p) < Kn(|« — v\). Property (h) follows 

from a straightforward analysis of the differential equation (e) which implies that 
^m(O) = for to > i. In the next section, we will see that the cusp at x — and 
the convexity of V m {x) in x > return when V m {x) is replaced by V^{x) as in the 
Slater model. 

The convexity of 1/V m (a;) can be rewritten as 

1 1 1 

\Vrn W) ~ V m (x) + Vjy) 

Using property (g) with a — ^, one easily finds that the this implies 

f 1 f 
< 1 . 

V m (x + y) V m (x) V m (y) 

Since 1/V m (x) « |a;| for large |a;|, this subadditivity inequality plays the role of the 
triangle inequality in applications. The proof of (i) is extremely delicate. Because 
1/V m (x) « x for large x, we need to prove the convexity of a function that is nearly 
linear so that its second derivative is extremely close to zero. Proving that this 
derivative is positive is equivalent to proving some rather sharp inequalities on the 
ratio V m (x)/V m -i(x). 

In the special case m = these inequalities (which are also discussed in [ BR1| 
and |SW|) are equivalent to 



(4.1) 

for x > 0, where 
(4.2) 



g*(x) < V (x) < g^x) 



9k(x) 



k 



(k - l)x + yjx 2 + k 



Multiplying (4.1) by x 



1 



obtain general ratio bounds, define 



converts this to a bound on the ratio 



Vo(x) 

v-i(xy 



To 



(4.3) 



ky 



(k — l)y — to + y/(y + to) 2 + ky 
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Then it is shown in RW| that 
(4.4) 



G%-\x 2 ) < J m{X ) , <Gf(x 2 ) 



V m -i(x) 

for all integer to > and x > 0. The sense in which these bounds are optimal is 
discussed in |RW], Our proof of these inequalities relies on an inductive argument 
and, hence, is valid only for integer to. A proof extending them to general real 
to > — 1 would immediately imply that properties (i) and (j) also hold for general 
real m > — 1. 

Another interesting open question is whether or not V m (x) is convex in m? In 
particular, is 2V m (x) < V m +i(x) + V m -\(x)l It follows from property (e) that this 
is equivalent to asking if V^ n (x) is increasing in m. 

5. Recursion and Averaged Potentials 



Using integration by parts on (1.3) one easily finds that V m satisfies the recur- 
sion relation 



(5.1) 



1 



V m (x) = — [(to - i - x 2 )U m _i(x) + x 2 V m - 2 (x)] 



for all to £ R, to > 1. Iterating this, one finds that when to is a positive integer 



(5.2) 



V m (x) = 



1 

2m 



m-2 



(1 - 2x 2 )V m - 1 {x) + \ V k {x) + 2\x\ 



k=0 



These relations are useful for studying V^(x). For example, it follows imme- 
diately from (5.2) that 



(5.3) 



2x 2 

Kv (x) - 2V N {x) - — \V-x{x) - V N -i(i 



This can then be used to show that Vj^(x) is convex for all x > 0. Furthermore 

lim — VJ^tx) = — -^r, verifying that VJ^(x) has a cusp at x = 0. 
x^o+ dx N 

It is interesting to note that (5.2) also implies that there are polynomials P m (y) 
and Q m (y) of degree to such that for integer to > 1 

(5.4) V m {x) = P m {x 2 )V {x) + xQ 



These polynomials have many interesting properties. In [RW] it is shown that 

1 



(5.5) 



P m (y) 



iFi 



,2' 2 



771,2/ I 



to B(m, i) 

where B(m,n) denotes the beta function and ii 7 i(a,7,?/) denotes the indicated 
confluent hypergeometric function. 
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